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Abstract 

We develop an extension of chiral perturbation theory to the A(1232)-resonance energy region 
and apply it to investigate the pion electroproduction off the nucleon (e“ N —> e~ N tt). We 
present a complete calculation of this process, in the A-resonance region, up to next-to-leading 
order in the 5-expansion. At this order, the only free parameters are the three low-energy constants 
corresponding to the magnetic (Ml), electric (E2), and Coulomb (C2) yN A transition strength. 
After fitting these parameters to a few well-known data, our calculation provides a prediction for 
observables and multipole amplitudes of pion electroproduction. These results compare favorably 
with the phenomenological multipole solutions and recent experimental results from MIT-Bates 
and MAMI. Our prediction for the pion-mass dependence of the yNA form factors offers an 
explanation for the discrepancy between the recent lattice-QCD results and the experimental value 
for the “C2/M1 ratio” at low Q^. 
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I. INTRODUCTION 


The first excited state of the nucleon — the A(1232) resonance — dominates pion- 
production phenomena and plays an important role in our understanding of the low-energy 
nucleon structure. High-precision measurements of the nucleon-to-A transition by means of 
electromagnetic probes became possible with the advent of the new generation of electron 
scattering facilities, such as BATES, MAMI, and JLab, many measurements being completed 
in recent years [ll 12 , y, 1 ^ • 

The electromagnetic nucleon-to-A (or, in short 7 AA) transition is predominantly of the 
magnetic dipole (Ml) type. In a simple quark-model picture, this Ml transition is described 
by a spin flip of a quark in the s-wave state. Any d-wave admixture in the nucleon or the 
A wave-functions allows for the electric (E 2 ) and Coulomb ((72) quadrupole transitions. 
Therefore by measuring these one is able to assess the presence of the d-wave components 
and hence quantify to which extent the nucleon or the A wave-function deviates from the 
spherical shape, i.e., to which extent they are “deformed” [^. The d-wave component 
of A’s wave-function can be separately assessed by measuring its electric quadrupole mo¬ 
ment. However, this would be extremely difficult because of the tiny lifetime of the A. 
The 'jNA transition, on the other hand, was accurately measured in the pion photo- and 
electro-production reactions in the A-resonance energy region. The E2 and C2 transi¬ 
tions were found to be relatively small at moderate momentum-transfers the ratios 

Rem = E2/M1 and Rsm = (72/Ml are at the level of a few percent. 

Traditionally , t he resonance parameters are extracted by using unitary isobar mod¬ 
els H Q i i Q □ in , which in essence are unitarized tree-level calculations based 
on phenomenological Lagrangians. However, at low the yAA-transition shows great 
sensitivity to the “pion cloud” jwhich until recently could only be comprehensively studied 
within dynamical models 0, Q llllii, E 0- (Unlike the isobar models, dynamical 
models include quantum effects due to pion loops.) 

With the advent of the chiral effective field theory (yEFT) of QCD 0,0 and its 
extensions to the A-resonance region 0,01, it has become possible to study the nucleon and 
A-resonance properties in a profoundly different way. Recently, we have been able to perform 
first yEFT studies of the yAA-transition in pion electroproduction 0] and of the A- 
resonance magnetic moment in the radiative pion photoproduction 0. The advantages over 
the previous dynamical approaches are apparent: yEFT is a low-energy effective field theory 
of QCD and as such it provides a hrm theoretical foundation, with all the relevant symmetries 
and scales of QCD built in consistently. Moreover, we hnd that already at next-to-leading 
order (NLO) in the “^-expansion” of Ref. 0 , the observables for pion electroproduction 
in the A-resonance region are described remarkably well. The yEFT, therefore, provides a 
theoretically consistent and phenomenologically viable framework, which, in particular, will 
allow for a model-independent extraction of the resonance parameters. 

Tremendous progress has recently been achieved as well in the lattice QCD simulations 
of the yAA transition. The present state-of-the-art results 0] are “quenched” and are ob¬ 
tained for pion masses above 300 MeV. These results can only be confronted with experiment 
after an extrapolation down to the physical pion mass of 140 MeV. A linear in the quark 
mass (rUg ~ m^) extrapolation was used in Ref. 0 . The thus obtained Rsm ratio, at low 
Q^, was found to be in major disagreement with experiment. The apparent caveat of this 
result is that the extrapolation in the quark mass needs not to be linear. The non-analytic 
dependencies, such as ArrH and InrUg, are known to be important as one approaches the 
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small physical value of niq. These non-analytic terms can also be obtained from yEFT, 
and, as we have demonstrated 0 , the Rem and Rsm ratios do exhibit a pronounced non- 
analytic quark-mass dependence, such that the lattice results 0 
experiment. Here we shall rehne our calculation by consistently including the 
dependence of the nucleon and A-isobar masses in the same yEFT framework 


can be reconciled with 
nark-mass 


A brief account of this work has recently been published j23[. In this paper we present 
an extensive description of the work, improve on our theoretical error analysis, and present 
some new results that could not be included in the brief publication. This paper is organized 
as follows. 

In Sec. ini we recall the relevant chiral Lagrangians, while Sec. IIIII explains the power 
counting in the 5-expansion scheme based on which the leading- and next-to-leading-order 
contributions are selected. The chiral-loop contributions to the yAA-transition form factors 
are evaluated in Sec.|TV|by using two different techniques. In Sec.|V|we discuss the theoretical 
uncertainty of our calculation due to the neglect of higher-order effects. In Sec. EH we 
present the results for pion photo- and electroproduction observables, multipoles, and the 
extracted yAA form factors. In that section we also discuss the yEFT predictions for 
the m^-dependence of the yAA transition and compare them with available lattice results. 
Sec. EH lists the main points and conclusions of the paper. The two Appendices contain 
technical details concerning the Feynman-parameter and dispersion integrals, respectively. 


II. THE EFFECTIVE CHIRAL LAGRANGIAN 

In this section we dehne the relevant effective Lagrangian of low-energy QCD, where we 
restrict ourselves to the two flavor, isospin symmetric case (m„ = rrid = r rin). The guiding 
principle is the symmetry under the SU(2)2,xSU(2)/^ chiral rotations, see [l9|, l20|, |27| . The 
lowest-order Lagrangian of the Goldstone-boson (pion) isovector held 7r“ (a = 1,3) 

£?> = \Sl Tr[ + (t; + C/t) X1 = - ImJ+ 0(1“) ], (1) 

where 

U{x) = = u\x), with 7rV“ = ^ , (2) 

and the linear in U, explicit symmetry-breaking term is proportional to the quark (or, pion) 
masses: y = 2Bmq = -|- O(m^). Furthermore, at this lowest order, ~ 92.4 MeV is the 

pion decay constant and B is related to the scalar quark condensate as B = — (qq) /f^. In 
the notation 12 ^)^ the superscript stands for the number of derivatives of Goldstone-boson 
helds and insertions of their mass. 

The chirally symmetric Lagrangian for the nucleon isodoublet held, A = [p, n)^, can 
conveniently be written by using the SU(2) vector and axial-vector currents, 

I r“u“(a;) = ^. {u + u^d^u) = ^ r“7r^ + 0(7r^), (3a) 

ir“a“(a:) = {v) d^u - ud^u^) = (3b) 

and the chiral covariant derivative acting on the nucleon held, 

Zl^A = a^A+izrX^- (4) 
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In terms of these definitions the nucleon Lagrangian begins at 


( 5 ) 


where gA — 1.267 is the nucleon axial-coupling constant. 

In the presence of the electromagnetic held {A^, the charge of the pions is accounted for 
by making the “minimal substitution”: —>■ in the above expressions. 

Similarly, the proton charge is included by the minimal substitution in the chiral derivative 
as: Dfj^N —>• 77^77 -|- ie^(l -|- P)A^N. 

The Lagrangian for the spin-3/2 isospin-3/2 A-isobar can be written in terms of the 
Rarita-Schwinger (vector-spinor) isoquartet field^, = (A++, A+, A°, A”)"^, by using the 

vector and axial-vector currents Eq. © in the isospin-3/2 representation of SU(2). The 
corresponding generators T“ have the following matrix representation: 


= - 
3 


/ 0 73/2 0 0 \ 

73/2 0 1 0 

0 1 0 73/2 

Vo 0 73/2 0 / 

/ 0 -73/2 0 0 


= 


2i 


\ 


73/2 0 -1 0 

0 1 0 -73/2 

VO 0 73/2 0 

= diag(l, l,-7-l), 

and satisfy = 5/3. The chiral derivative is then given by: 

D^A, = d»^A, + *7r“A,, 

and the hrst-order chiral Lagrangian can be written as: 


4^ = a^ 


ii7AA^7r75A7 


( 6 a) 


( 6 b) 

( 6 c) 


(7) 

( 8 ) 


where ~ 1.232 GeV is the mass of the A-isobar, Ha is the axial coupling constant given, 
in the large-limit j^, by Ha = (9/5)qa and this value is known to be consistent with 
the empirical information, see, e.g., [29|. The totally-antisymmetric products of y-matrices 
are dehned as: 7 ^^ = The electric charge of the A is accounted 

for by the following minimal substitution: + ^e^l + dT^)A^. 

The free Rarita-Schwinger held obeys the following held equation: 


= Ma 7 ^"A„ (9) 

which, in particular, yields the constraints: 7 ^A^ = 0 = d^A^. The constraints reduce 
the number of spin degrees of freedom of the vector-spinor held to the physical number 
appropriate for a massive particle with spin 3/2. These constraints arise essentially due to 


^ Frequently chiral Lagrangians for the A-isobar are written in terms of the isovector-isodoublet field |2ll: 
A“ = T“A^, with T isospin-l/2-to-3/2 transition matrices. It is an alternative to the isoquartet represen¬ 
tation adopted in this work. The two representations, however, are equivalent at the level of observables. 
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the invariance of the free-held kinetic term [hh.s. of Eq. ©], upon the gauge transformation 
of the spin-3/2 held: 

Af,{x) ^ Af,{x) + df,e{x), ( 10 ) 


where e is a spinor. In order for interactions to support the number of the free-held con¬ 
straints, they must be symmetric with respect to the same (or, at least, similar) transfor¬ 
mation 0. Now, the chiral interactions in Eq. (jH)) obviously do not have this symmetry. 
Therefore, we redehne them by using the free held equation, or, equivalently, redehne the 


spin 3/2 held, such that they become symmetric 
coupling becomes 


[31|. Then, e.g., the lowest-order axial 


+ o(^=). (ii) 

2MA/7r 

The diherence with the original interaction term from Eq. (jH)) is of higher order (in the pion 
derivatives and mass). 

We write the iVA-transition Lagrangian right away in such an expanded form which 
exhibits the spin-3/2 gauge symmetry: 


C 

C 

c 


( 1 ) 

NA 

( 2 ) 

NA 

(3) 

NA 


iVT“ 7 ^"^ (a^A,) aA7r“ + H.c., 


2UMa 




2M]\f[M]^ -|- T/a) 
-3e 

2M]^[M]\f M/\ 




-ATS 


gEid^A^) -F ^7“(aoA^ - d^A^) 
Ma 


( 12 a) 

( 12 b) 

F'"' + H.c.,(12c) 


where and are the electromagnetic held strength and its dual, T“ are the isospin- 
l/2-to-3/2 transition (2 x 4) matrices: 


rjn2 


S (0 1 0 \ 

S V 0 -1 0 v/3S 

S Vs 0 1 0 \ 

S V 0 1 0 Vss 

/2 /o 1 0 0\ 

V 3 VO 0 1 OS 


(13a) 

(13b) 

(13c) 


satisfying T“T ^1 = 

Note that the electric and the Coulomb 7 AA couplings are of one order higher than the 
magnetic one, because of the 75 which involves the “small components” of the fermion helds 
and thus introduces an extra power of the 3-momentum. 

For the momentum-space A propagator we use 


5“^(p) 


|i -F Ma 

Mi -p 2 


0/3 


+ 


(i-C)(Ci^ + ^a) . a 3 

3{CV-Ml) 


7V) + 


2(1 -C^)pV’ 

3{CV-Ml) _ 


1 


where C is a gauge-hxing parameter, see 
gauge for spin -1 case here is, V = 00 : 


(14) 

Refs. 00 for details. The analog of the Landau- 




f + Ma 
Ml-p^ 


p{3/2)o/3^p)^ 


(15) 
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with 


-p(3/2) 



p2 J 


3p2 


Pm, 


(16) 


the covariant spin-3/2 projection operator. As long as the conpling of the spin-3/2 held are 
gauge-symmetric with respect to the transformation (HDD, the results are independent of the 
choice of the gauge-hxing parameter (. 


III. POWER COUNTING AND RENORMALIZATION 


The inclusion of the A-resonance introduces another light scale in the theory, the reso¬ 
nance excitation energy: A = Ma — Mn ~ 0.3 GeV. This energy scale is still relatively light 
in comparison to the chiral symmetry breaking scale ~ 1 GeV. Therefore, 6 = A/A^^sb 
can be treated as a small parameter. The question is, how to compare this parameter with 
the small parameter of chiral perturbation theory (yPT), e = m.,^/A^ sb- After all, our aim 
is to organize an expansion in a small parameter and to estimate the size of various contri¬ 
butions based on power-counting rules. So some specihc relation between 5 and e would be 
helpful. 

In most of the literature (see, e.g., 0 ) they are assumed to be of comparable size, 5 ~ e. 
This, however, leads to a somewhat unsatisfactory result: the A-resonance contributions 
are always estimated to be of the same size as the nucleon contributions. Thus, obviously, 
the A-contributions are overestimated at lower energies and underestimated at the resonance 
energies. To estimate the A-resonance contributions correctly, and depending on the energy 
region, one needs to count 5 and e differently. 

A relation e = was suggested and explored in Ref. 0, and is adopted in this work. 
The second power is indeed the closest integer power for the relation of these parameters in 
the real world. We should stress that this relation is used for power-counting purposes only. 
It is not imposed in the actual evaluations of diagrams. Each diagram is simply characterized 
by an overall ^-counting index n, which tells us that its contribution begins at 0(5"’). 

Because of the distinction of and A the counting of a given diagram depends on 
whether the characteristic momentum p is in the low-energy region [p ~ mjr) or in the 
resonance region (p ~ A). In the low-energy region the index of a graph with L loops, 
pion propagators, nucleon propagators, Aa A-propagators, and V vertices of dimension 
i is 

2 ( + 4L - Aa, - 2AJ - Aa = 2n^pr - Aa, (17) 


n = 


where n^p-r is the index in xPT with no A’s j27|. In the resonance region, one distinguishes 
the one-A-reducible (OAR) graphs see e.g., graph (a) in Fig. |21 Such graphs contain 
A propagators which go as 1/ (p — A), and hence for p ~ A they are large and all need 
to be included. This gives an incentive, within the power-counting scheme, to resum A 
contributions. Their resummation amounts to dressing the A propagators so that they 
behave as l/(p — A — A). The self-energy A begins at order p^ and thus a dressed OAR 
propagator counts as 1/5^. If the number of such propagators in a graph is Nqar, the 
power-counting index of this graph in the resonance region is given by 


n — n^pT — Aa — ^Nqar, 
where Aa is the total number of A-propagators. 


( 18 ) 
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A word on the renormalization program, as it is an indivisible part of power counting in a 
relativistic theory. Indeed, without some kind of renormalization the loop graphs diverge as 
A-^, where A is an ultraviolet cutoff, and A/” is a positive power proportional to the power¬ 
counting index of the graph. Also, contributions of heavy scales, such as baryon masses, 
may appear as M-^. The renormalization of the loop graphs can and should be performed so 
as to absorb these large contributions into the available low-energy constants, thus bringing 
the result in accordance with power counting 0. 

To give an example, consider the one-7riV-loop contribution to the nucleon mass. Fig. 
For the ttNN vertex from the power counting tells us that this contribution begins 
at O(m^). An explicit calculation, however, will show (e.g., |0) that the loop produces 
0(m°) and O(m^) terms, both of which are (inhnitely) large. This is not a violation of 
power counting, because there are two low-energy constants: the nucleon mass in the chiral 
limit, and Cim, which enter at order 0(m°) and 0{rnff), respectively, and renormalize 

away the large contributions coming from the loop. The renormalized relativistic result, up 
to and including O(m^), can be written as |0: 


Mn = 


- 4 ciN ml 
^9 a 


(8vr/^ 


■m^ 


4 1- 


+ 


AMI 


mn 

8Mj^ 


5/2 


arccos 


30 - 10 


m^ 


M. 


N 


2M]\f 

2 

+ 


+ 




17m^ 

lOAf/v 

4' 


M, 


N 




2 M, 


N 


In 


m^ 


M, 


N 


(19) 


and one can easily verify that the loop contribution begins at O(m^) in agreement with 
power counting. 

We now turn to the analysis of the pion electroproduction process. The pion electropro¬ 
duction amplitude to NLO in the 5 expansion, in the resonance region, is given by graphs 
in Fig. a) and (b), where the shaded blobs in graph (a) include corrections depicted in 
Fig.|2fc-f). The hadronic part of graph (a) begins at 0{6^) which here is the leading order. 
The Born graphs Fig. m contribute at 0{6). We note that at NLO there are also vertex 
corrections of the type (e) and (f) with nucleon propagators in the loop replaced by the A- 
propagators. However, after the appropriate renormalizations and <C A A, these graphs 
start to contribute at next-next-to-leading order. 

The A self-energy Fig. |2^c) can, in the ( = oo gauge, be presented as 




( 20 ) 


where has the spin-1/2 Lorentz form, computed by us earlier Hi Hi- To recall these 

results, it is convenient to introduce the following dimensionless quantities: 


Mn . A 

= 

(3 = \{l - r"^ + p?), a = 1-/3, (21) 

A2 = 1(52 _ + rf - /i2] = (3^-p^ = a^- r2. 


The imaginary part of the A self-energy is related to the resonance decay width, which at 
the leading order (LO) gives: 

Fa = -2 ImrA(MA) = {An/?>) M^ClX\a + r). (22) 
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with Ca = hyiMA/(87r/jr). The experimental value for the A-resonance width, Fa = 
0.115 GeV, hxes Ha — 2.85. At NLO, the residue of the A propagator receives a correction: 


Im A^(Ma) = —2ttC\ X \a (3 {a-\- r) — gA^(r + — /i^)] . (23) 

The renormalized A-propagator, in the (^ = oo gauge, at NLO reads: 

^ - Ma)[ 1 - * Im r^(MA)] + \ir^ • 

The real part of the renormalized A self-energy, Ta, contributes to the mass as: 

Ma = Mf - 4 ciA ml + (Ma), (25) 

where 


^a(Ma) = -ICI Ma [V{fi, 6) - G(0, 6) - ^(o, 6)] 

The ttN loop integral V (p, 5) is given in analytical form by jlf 


(26) 


V (p, S) = ■|(r -f a) [/3 (/i^ — 2A^) Inp^ -|- a (r^ 
+ 4A2 + 4A^ 0(A) ] + A/(ln/i2 - 1) - 


with the elementary function O defined as: 


0(A) 


_L 1 ^ 

2X 13-fi'^+X ’ 


1 


arctan 


a/3+A2 


- 2A^) Inr^ - |(q;^ -f (3^) 


. ir‘^(lnr^ - 1), 

(27) 

A^ > 0 


. A^ < 0 

(28) 


The region A^ > 0, corresponds with < Ma — M^r, where the A is unstable. The region 
A^ < 0, corresponds with tiIt, > Ma — M^, where the A is stable. 

Furthermore, from Eq. (EH) we read off the ml and ml terms, which enter Eq. (EHD and 
need to be absorbed by the renormalization of the low-energy constants. 


1 /( 0 , 5) 


E'(0,5) 


2 " 

^( 2 - 0 + 2(1 


Inr^ — :j^(l + J’)^(l — lii(l — r^) 



(29) 


(30) 


The vector-meson diagram, Fig. Hd), contributes to NLO for ~ AA. We include it 
effectively by giving the (^M-term a dipole Q^-dependence (in analogy to how it is usually 
done for the nucleon isovector form factor): 


9m 


__ 

(l + gVO-71 GeVy 


(31) 


Analogous effect for the qe and gc couplings begins at N^LO. The loop corrections to the 
yAA-vertex are discussed separately, in the following section. 















IV. EVALUATION OF THE -/NA FORM FACTORS 


In this section we present a detailed analysis of the yVA-vertex corrections that appear 
at NLO, see Fig. Efe ) and (f). To the order we work, the pseudovector irNN interaction 
from the Lagrangian © is equivalent to the pseudoscalar one, and we use the latter in the 
actual calculations of the graphs Fig. Etc, f). 

The rest of the section is organized as follows. First we will comment on the electro¬ 
magnetic gauge invariance of these contributions. Then, after a brief summary of different 
decompositions of the yiVA vertex into Lorentz-invariant form factors, we evaluate the chiral 
corrections to those form factors at NLO. As a cross-check, we do the calculations using two 
different techniques: the Feynman-parameter method and the sideways dispersion relations. 


A. Electromagnetic current conservation 

The addition of the one-loop graphs Fig. He, t) to the “tree-level” 'jNA vertex from 
Eq. ()12|1 . should not spoil the e.m. current conservation: 

Qt, Ua{p) «(p) = 0 ) (32) 

where Ua is the free Rarita-Schwinger vector-spinor of the A, u is the free Dirac spinor of the 
nucleon, and momenta are defined as in Fig. El This sort of conditions are most conveniently 
checked by using the Ward-Takahashi (WT) identities for the e.m. couplings of the nucleon, 
pion and A fields: 


Qv = dNq-l = eN [S^^ip) - S/{p )], (33a) 

Qu = e^q-{p + p') = (p') - (p)], (33b) 

qu = eA = eA[^Z'“^(p') - 5 Z'“^(p)] , (33c) 

where q = p' — p, cn = + T 3 ), e^^ = iee°'Yc, ca = ^e(l -|- 3T^) and 5'“^ denote the 

corresponding inverse propagators. Applying these identities to the diagrams Fig. He. t) we 
obtain: 

qi^Uaip) F“)(f^M(p) = eM„(p') [A^a(p) “ ^naY)] Yp), (34) 

where A^a one-loop A —>■ A self-energy, which, due to Lorentz covariance has the 

following general form: 

s^Yp) = Mp")p^ + b{p^)Y- (35) 

The scalar functions A and B do not vanish, however we may use the Rarita-Schwinger 
conditions, p ■ u{p) = 0 = 'j ■ u{p), to show that the second term on the r.h.s. of Eq. 
vanishes. The hrst term can only be canceled by the contribution of the diagram in Fig. EJ 
Thus, the condition of the e.m. current conservation requires the diagram of Fig. E] to be 
included in addition to the NLO contributions Fig. Efe, f). When computing this diagram, 
we discover that it is actually of N^LO, because the Lorentz form Eq. (ESD cuts out the spin- 
3/2 part of the A-propagator and hence also the pole of the propagator. The A-propagator 
in Fig. m therefore, counts as 5°. 

To NLO, the sum of Fig. He) and Fig. m) conserves the e.m. current. 
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B. Form-factor decompositions 


The 'jNA vertex can in general be decomposed into three Lorentz covariants. For in¬ 
stance, 


Ua{p) 


X 


1^3 ii^A T dT/v 

V 2 M^[(Ma + Mjv )2 + Q 2 ] 


Ua{p') { gAiiQ^) p'^ qx 

+ 9e{Q^) {q'^P^ -q- P *75 
+ £/c(Q^) £/"'") * 75 } u{p), 


(36) 


where q = p' — p is the photon 4-momentnm, Q‘^ = —q^. Furthermore, gM^Q"^), gE^Q^), 
and gc{Q^) are the magnetic dipole, electric quadrupole, and Coulomb quadrupole form 
factors, as defined in Ref. [l^. In the limit = 0 they are equal to the physical values 
of the corresponding parameters in the Lagrangian m- These form factors relate to the 
conventional magnetic (G^), electric (G^) and Coulomb (G^) form factors of Jones and 
Scadron as follows: 


G*m = gu + {-P'y gE + Q^gc) , 

v+ 

^*E = -qi-gE + Q^gc), (37) 

2 M\ 

Gc — — yyr' (Se + P7 fi'c), 
v+ 


where Q± = ■>/ (Ma ± = Q^/^ai /^7 = ^{7 — — Q'^)■ The multipole ratios 

E2/M1 and G2/M1 at the A-resonance position can be expressed in terms of these form 
factors as: 


Rem — 


G 1 


Rsm = 


Q+Q- G 


c 


G*j^ 4M_^ G%i 

Alternatively, the on-shell 7 AA vertex is also often expressed in the following form : 


(38) 


u„(p)T'^%^u(p) = ^Jlua(p) {( 7 '^g" - 7-?£/“^) ^i(Q^) 

+ (q-p'g^^-q'^p'ng2(Q^) (39) 

+ (g“ gz(Q‘^) ]iy^u(p), 


which is defined for the p —> A"^ transition, yielding the isospin factor \f 2 J 3 . The form 
factors gMiPE and gc can then be expressed in terms of gi(Q^)(* = 1,2,3) through the 
relations : 


gniQ"^) 

gE(Q^) 

gc(Q^) 


MmQI 


3 M^(M/x T TJ/v 
2 MmQI 


■gi(Q^ 


3 + M}^ 

2 


{9i(Q2)-MaJ2(7)}. 


3 (ATa T Mj\f) 


g^iQ')- 


(40) 
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C. Feynman-parameter method 


The one-loop corrections to the jNA form factors are given by the graphs in Fig. He) 
and (f). As they contain three propagators, we apply the Feynman-parameter trick : 

[At, + B{x - y) + C{1 - x)]^' 

For the diagram Fig. H®) we have A = + ie, B = {k + + ie^ C = 

(p — -t-ie, where k is the integration 4-momentum. After the shift of the integration 

momentum, k —>■ k — {x — y)q -l- (1 — x)p, we obtain : Ay + B{x — y) + C{1 — x) —>• 
where 


M‘^{x, y) = im?^x+ (Af^ — p^a;)(l — x) — 2p ■ q{l — x){x — y) — q^{x — y){l — x — y) — ie. (42) 

For the case of the N ^ A transition we may use: = M^, p'^ = (p + q)'^ = and 

hence 2p ■ q = M\ — and 

M^{x, y) = mix + [M% - xM\ + {M\ - M^)p] (1 - a:) -h q^y{x - y) - ie. (43) 


The next step is to perform the integration over 4-momentum which can be done by using 
the following rules: 


' 1 _ , , 2^ .(-l)”F(n-2) 

(27r)4 (P - m?Y ~ ^ * (47r)2 F(n) "" 


d^k 


k^ky 


(27r)^ {}A‘ — m^)” 2{n — 1) 


) Qfiv ) 


(44a) 

(44b) 


Integrals with an odd number of 4-vectors k in the numerator vanish. The integral Jn 
diverges for n = 1 and 2, but can be dehned via dimensional regularization as: 


Ji(m^) 

J2{m^) 


—im 

(47r)2 

—i 

(47r)2 


2 

2 , 




(45a) 

(45b) 


where d —> 4“ is the number of dimensions and je = ~r'(l) ~ 0.5772 is the Euler constant. 

Then, the (MS'-subtracted) result for the graph of Fig. He), where the photon couples to 
the charged pion in the loop, can be decomposed into the three Lorentz covariants Eq. dnni) 
and cast in the form : 


1 i-y 

= -Cna J dyy j dx InAl^, (46a) 

0 0 

1 i-y 

= +Cna j dyy J dx {\nM‘l - 2x [(1 - p) (1 -t- r) - x] M~‘^] , (46b) 

0 0 

1 i-y 

= -Cm A j dyy{2y - 1) j dx [(1 - p) (1 + r) - a:] (46c) 

0 0 

with Ail = [x — PY — -|- 2f3'yxy + QYyi)- — y) — is, (46d) 
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while the analogous contribution of the graph Fig. m), where the photon couples to the 
electric charge of the nucleon in the loop, is given by : 


^-y 


9m = -Cna j dyy j dx {2lnMj 
0 0 

+ Mj'^ [2a^xy + Q^y{l - y) - {x + r){l - x)] } , 


1 i-y 


9e^ = -Qm +2Cna I dyy I dxM/ [{x + r){l-x)-xy{l+r)], 


M) 


= -2Cna I dyy^ I dx[y{l + r) - 1 + x]Mf'\ 


0 0 
i-y 

,2 


0 

with A4‘j = (x — a)^ — + 2ajxy + Q^y(l — y) — is, 

where we have used the dehnitions of Eq. (EH) and introduced, 

QW 


Cna = 


(87r/^)2 3(l + r)' 


(47a) 

(47b) 

(47c) 

(47d) 

(48) 


Finally, the integration over x and some of the r/-integration can be done analytically, see 
Appendix A. 


D. Dispersion method 

Alternatively, we can compute the correction to the yiVA vertex by exploiting the ana- 

lyticity of the loop contributions. We use the so-called sideways dispersion relations jH, , 

and hence start with calculating the absorptive (imaginary) part of the yiVA vertex from 

the cut in the 7rA-loop diagrams of Fig. El Subsequently, the real part of the on-shell yAA 

vertex is computed through the dispersion relation in the A-resonance 4-momentum squared 
/2 

p . 

We hrst express the half off-shell yAA vertex (where we allow the A to be off-shell) in 
terms of the invariants gi{s = Q^), for i = 1, 2, 3, as : 

f“^m(p) = (^py^(p)) ^ { (tV - gtis, Q^) 

+ {q-pg^'"-q^p^) gt{s, Q^) 

+ [q^ qf^ - q'^ g^^^) g+{s, Q^) }i^5u{p) 

+ (p5^(p)) ^ { (aV - i-qg^^) 9i{s, Q^) 

+ {q-p g^^ -q'^p^) g 2 {s, Q^) 

+ (/ q^ - q^ /'") 9z ^ 75 u{p), (49) 

where the spin-3/2 positive (negative) energy projectors are given by 
(^z/2f ,A (y-p'±MA) [ 1 I ^ 

=- 2 Ml- 1 ®“'’ ■ 3'°'’^ “ 3 ^ ® ® 7 ' * ' 
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We next determine the real parts of the invariants for s = M\ (i.e., the A is on-shell) 
through a dispersion relation in s as : 


Regt{Ml,Q^) 




Sth 


Im^+(g,Q^) 

s-Ml ’ 


(51) 


where the integration starts from the vriV threshold Sth = (d^Af + ’^ 7 r )^5 and where P denotes 
the principal value integration. These unsubstracted dispersion integrals do not converge. 
However, the dispersion relations for the following differences 


Re^+(Mi,g2)_Re^+(Mi,0) 




^th 


s-Ml 


(52) 


do converge, for all the considered contributions. This subtraction can be put in cor¬ 
respondence with the renormalization of appropriate low-energy constants from the La- 
grangian m- 

The expressions for the absorptive (imaginary) parts of these form factors are calculated 
in Appendix IbI with the results given by Eqs. (IB50f) . ()B51jl . The dispersion integrals in 
Eq. are evaluated numerically. The physical form factors are then simply given by 
g,m=9tiMl,Q^). 

We compared the thus obtained results with the Feynman-parameter method results. We 
found perfect agreement between both methods for all three jNA transition form factors, 
which provides a cross-check on our calculations. 


V. ERRORS DUE TO NEGLECT OF N^LO EFFECTS 


Prior to presenting the results of our NLO calculation, we would like to make an estimate 
of the theoretical uncertainty due to the neglect of higher-order effects. Some of the next- 
next-to-leading order (N^LO) contributions are shown in Fig. El Of course, there is no 
substitute for an actual calculation of those effects, but at present we only know that they 
must be suppressed by at least one power of S (= A/A^sb) as compared to the NLO and two 
powers of 6 as compared to the LO contributions. Therefore, we can estimate the size of the 
N^LO contribution to an amplitude A as: Anlo or Alo The theoretical uncertainty of 
a calculation up to and including NLO can thus be estimated as: 

Aerr = |^LO + ^NLO I ■ (53a) 


In cases where the amplitude does not receive any LO contributions, we have no other option 
than, 

Aerr = |AnLo| (53b) 

This looks nice and simple, however, there are a few caveats in the implementation of 
such an estimate. First of all, although we have introduced 6 as A/A^sb, h counts as well 
the other small scales of the theory, tUtt and Q^. Therefore, we shall estimate the error using 
the following expansion parameter (assuming A^sb ~ Mjsf): 


6 


1 

3 


A 



Mj,) 


(54) 
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where all the light scales are treated on equal footing and hence are averaged over. 

Secondly, what if the amplitude happens to vanish at some kinematical point. According 
to Eq. the theoretical calculation at that point would be perfect, which is of course 
unlikely to be true in reality. So, when considering dependencies on kinematical variable(s), 
we shall take an average of the error over some appropriate region of that variable. 

Given these two points, we are led to the following formula for the theoretical uncertainty 
of the NLO calculation for an amplitude A, 

\ |AUJ, LO = 0, 



and the subscript “ah’ indicates that the appropriate averaging is performed. 
The theoretical uncertainty of the NLO calculation of an observable O is: 


r 2\0\a,~S^ LO^O 

\ 2\0\aJ, LO = 0, 


(56) 


where the factor of 2 takes into account that an observable is a product of two amplitudes. 

Note that this error estimate differs from the one presented in our earlier paper j^, 
where we overlooked the case of vanishing LO contributions. This underestimated the error 
in some observables, as will be discussed in more detail below. 


VI. RESULTS AND DISCUSSION 

We are now in position to discuss the NLO results for pion photo- and electroproduction 
amplitudes and observables. We begin with the mult ip ole-analysis for pion photoproduction. 
The resonant photoproduction multipoles are well-established, and we use them to determine 
the two low-energy constants: Qm and the strength of the Ml and E2 yVA transitions. 
Subsequently, we discuss the results for pion e/ecfroproduction: cross sections, multipoles and 
the dependence of the Rem and Rsm ratios. The comparison with pion electroproduction 
observables allows to determine the third low energy constant in the xEFT framework, 
related to the strength of the C2 'jNA transition. Once the three low energy constants are 
determined, the Q^- and m 7 r-dependencies follow as a prediction of the NLO result. We 
discuss the predictions for the 771 ,^-dependence of the magnetic dipole jNA form factor, the 
Rem and Rsm ratios, and compare with the recent lattice-QCD results. 


A. Pion photoproduction 


In Fig. [71 we show the result of the yEFT calculations for the pion photoproduction 

f3/2) f3/2') 

resonant multipoles and Elf , around the resonance position, as function of the 

total c.m. energy W of the tiN system. These two multipoles are well established by the 


MAID 


nil and SAID [3^ partial-wave solutions which allow us to ht the two low-energy 
constants of the chiral Lagrangian Eq. m as : qm = 2.9, qe = —1.0. As is seen from the 
hgure, with these values the NLO results (solid lines) give a good description of the energy 
dependence of the resonant multipoles in a window of 100 MeV around the A-resonance 
position. Also, these values yield Rem = —2.3 %, in a nice agreement with experiment [ij. 
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The dashed curves in Fig.[7|show the contribution of the A-resonant diagram of Fig.|21^a) 
without the NLO vertex corrections Fig. 06, f). For the Mi+ multipole this is the LO 
contribution. For the i?i+ multipole the LO contribution is absent [recall that the qe coupling 
is of one order higher than qm, see Eq. m]- Hence, the dashed curve represents a partial 
NLO contribution to therein. 

Upon adding the non-resonant Born graphs Fig. Hb) to the dashed curves we obtain the 
dotted curves in Fig. 0 These non-resonant contributions are purely real at this order and 
do not affect the imaginary part of the multipoles. One sees that the resulting calculation 
is flawed because the real parts of the resonant multipoles now fail to cross zero at the 
resonance position and hence unitarity, in the sense of Watson’s theorem [^, is violated.^ 
The complete NLO calculation, shown by the solid curves in the hgure includes in addition 
the vertex corrections Fig. He, f), which restore unitarity exactly. This may come as a 
surprise, since we are doing a perturbative calculation, without a resummation of rescattering 
contributions. Nevertheless, it not difficult to see that our numerical result is not a fluke. 
Watson’s theorem is satisfied exactly by the NLO, up to-one-loop amplitude given the graphs 
in Fig. 0 

It is also interesting to examine our calculation for the non-resonant multipoles, which 
all receive contributions of NLO only. In Fig. 0 we show the NLO calculations for the real 
parts of the non-resonant s-, p- and d-wave pion photoproduction multipoles in the A(1232) 
region in comparison with the two state-of-the-art phenomenological multipole solutions, 
MAID and SAID. Note that at NLO in the h-expansion, the non-resonant multipoles are 
purely real. The multipole solutions show indeed that the imaginary parts of non-resonant 
multipoles, around the A resonance, are negligibly smaller than their real parts. From 
the hgure one thus sees that for most of the non-resonant multipoles, the parameter-free 
NLO results, agree fairly well with the phenomenological multipole solutions. The largest 
differences are observed for multipole. This multipole corresponds with nucleon 

quantum numbers. The cause of the appreciable difference in this channel is largely due 
to the nucleon anomalous-magnetic-moment contributions, which are not included in our 
calculation (since they appear at N^LO in the h-expansion), but which are included in the 
phenomenological solutions. 


B. Pion electroproduction 

We now will present the electroproduction observables obtained from the NLO amplitude 
of Fig. 0 The hve-fold pion electroproduction cross section can be expressed as : 

da da 


^ Recall that Watson’s theorem is a simple statement of unitarity based on a coupled-channel scattering 
equation to leading order in the electromagnetic interaction. The theorem relates the phase of a photo¬ 
production multipole to a corresponding (in spin and parity) 7riV-scattering phase-shift : 

For the resonant channel the corresponding phase-shift crosses 90 degrees at the resonance position, and 
thus the real part of the resonant multipoles must vanish. 


(57) 
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where the virtual photon flux factor is defined as : 


r E'^ (FF^ - Ml) 1 

" (27r)3Ee 2Mn (1-£)Q2’ ^ ^ 

where W is the invariant mass of the final tiN system, E^. {Ej) are the initial (final) electron 
lab energies, and £ denotes the photon polarization parameter. 

The 7 *A^ —> ttN cross section for unpolarized nucleons is expressed in terms of 5 response 
functions as : 


da 


dar ^ daL ^ . darr 

+ £ - h £ COS 2$ 


dVt^ 


dVt^ 


dVt^ 


\/2e{l + e) cos<h + h\/2e{l - e) sin<h 


(59) 


where and <h are the pion polar and azimuthal c.m. angles, respectively, and h denotes 
the electron helicity. 

In Fig.IHlwe show the NLO results for the different virtual photon absorption cross sections 
entering Eq. (El at the resonance position, and for ~ 0.127 GeV^, where recent precision 
data are available. Besides the low-energy constants qm and qe-, which were fixed from the 
resonant multipoles in Fig. H the only other low-energy constant from Eq. dH entering the 
NLO electroproduction calculation is go- hi Eq. (El, the main sensitivity on gc enters in 
aLT- A best description of the air data in Fig. IHlis obtained by choosing gc = —2.36. 

The theoretical uncertainty of the NLO result is estimated by using Eq. (El, where the 
average is taken over the range of ©jr- Note that air and aLT' do not receive any LO 
contributions and therefore the LO= 0 case in Eq. (El must be applied in the estimate. 
This point was overlooked in our first error estimates of the NLO calculation [^, which 
clearly led to an underestimate of the theory error for air and aiT'- 

From Fig. IHl one sees that the NLO xEFT calculation, within its accuracy, is consistent 
with the experimental data for these observables. 

The reliability of the present yEFT calculation can also be tested by comparing its 
predictions for the non-resonant multipoles with the phenomenological multipole solutions 
MAID and SAID. In Fig. CHI we show this comparison for the non-resonant s- and p- 
wave pion electroproduction multipoles at the resonance position as function of Q^. Note 
that there is a considerable uncertainty in the phenomenologically extracted s-wave scalar 
multipoles and at low Q^. As one can see from the figure, our calculation is in 

a reasonable agreement with the phenomenological solutions for most of the s- and p-wave 
non-resonant multipoles. The largest discrepancy is observed in the Q^-dependence of both 
and This discrepancy will need to be resolved by s-wave N^LO corrections 

which grow with Q^. 

In Fig. [TT]we show the dependence of the ratios Rem and Rsm- Having fixed the low 
energy constants gM, gs and gc, this dependence follows as a prediction. The theoretical 
uncertainty here (shown by the error bands) is estimated according to Eq. (I55|l . for LO=0 
case, and the average taken over the range of from 0 to 0.2 GeV^. From the figure one 
sees that the NLO calculations are consistent with the experimental data for both of the 
ratios. 

To see how higher-order effects may affect the dependence, we include the vector- 
meson type of dependence for the electric 'yNA transition by the following replacement of 
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the low-energy constant: 


9e 


( 60 ) 


(1 + QVhlf 

In contrast to the analogous effect for qm, which is of NLO, the inclusion of the de pen - 
dence in qe is a N^LO effect. We choose so as to satisfy the asymptotic condition |42| : 
Rem 1, for > cxo. This is achieved by the choice : O.TlGeV^. 

The dashed curve in Fisf. fTTl shows the resulting effect of the replacement (lt)()|l . as imposed 
on the complete NLO result shown by the solid curves with error bands. The fact that the 
dashed curves go out of the error bands at some point indicates that our error estimate is 
not designed for such high values. 


C. Chiral behavior and chiral extrapolations 

Since the low-energy constants qm, Qe, and gc have been hxed, our calculation can 
provide a prediction for the dependence of the yA^A transition form factors. The study 
of the m.n--dependence is crucial to connect to lattice QCD results, which at present can only 
be obtained for larger pion masses (typically > 300 MeV). 

The dependence of the nucleon and A-resonance masses, given above by Eq. (IT^ and 
(j20|l . are compared with lattice results in Fig.^J We constrain one of the two parameters in 
Eq. (jl0|l by the physical nucleon mass value at = 0.139 GeV, while the other parameter 
is £t to the lattice data shown in the hgure. This yields : = 0.883 GeV and Ci^v = 

—0.87 GeV“^. As is seen from the figure, with this two-parameter form for Mv, a good 
description of lattice results is obtained up to ~ 0.5 GeV^. 

Analogously to the nucleon case, we £x one parameter in Eq. (EHl) from the physical value 
of the A mass, while the second parameter is fit to the lattice data shown in Fig.^1 yielding : 

= 1.20 GeV and Cia = —0.40 GeV“^. As well as for the nucleon, this two-parameter 
form for Ma yields a fairly good description of the lattice results up to ~ 0.5 GeV^. 

In Fig. Uni we examine the m.n--dependence of the magnetic yVA-transition form factor 
in the convention of Jones and Scadron, see Eq. (jSZj). At the physical pion mass, this 

3/2 

form factor can be obtained from the imaginary part of the multipole at hF = Ma 

(where the real part is zero by Watson’s theorem) as : 

' ImMf/'(W = (61) 

V "J CK em IQ. I / 

where |p* | (|q* |) denote the pion (virtual photon) c.m. three-momenta respectively at the 
resonance position, i.e. for W = Ma. Recall that the value of G*^ at = 0 is determined 
by the low-energy constant gu- The Q^-dependence then follows as a prediction of the NLO 
result, and Fig. IT^ shows that this prediction is consistent with the experimental value at 

= 0.127 GeV^ and physical pion mass. 

The ruTT-dependence of is also completely hxed at NLO, no new parameters appear. 
In Fig. [131 the result for at Q‘^ = 0.127 GeV^ is shown both when the m^-dependence of 
the nucleon and A masses is included and when it is not. Accounting for the m^r-dependence 
in Mat and Ma, shown in Fig. apparently changes the result for quite significantly. 
The yEFT calculation, with the dependence of and Ma included, is in a qualitatively 
good agreement with the lattice data shown in the hgure. The yEFT result also follows an 
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approximately linear behavior in although it falls about 10 - 15 % below the lattice 
data. This is just within the uncertainty of the NLO results. One should also keep in 
mind that the present lattice simulations are not done in full QCD, but are “quenched”, so 
discrepancies are not unexpected. 

In Fig. HH we show the m^-dependence of the ratios Rem and Rsm and compare them 
to lattice QCD calculations. The recent state-of-the-art lattice calculations of Rem and 
Rsm 11^ use a linear, in the quark mass {nig oc m^), extrapolation to the physical point, 
thus assuming that the non-analytic m^-dependencies are negligible. The thus obtained 
value for Rsm at the physical value displays a large discrepancy with the experimental 
result, as seen in Fig. El Our calculation, on the other hand, shows that the non-analytic 
dependencies are not negligible. While at larger values of m^, where the A is stable, the 
ratios display a smooth dependence, at = A there is an inflection point, and for 
< A the non-analytic effects are crucial, as was also observed for the A-resonance 
magnetic moment mill. 

One also sees from Fig. HU that, unlike the result for G]^, there is only little difference 
between the yEFT calculations with the rriTr-dependence of Mat and Ma accounted for, and 
our earlier calculation j^, where the ratios were evaluated neglecting the m^-dependence of 
the masses. This is easily understood, as the main effect due to the dependence of Mat 
and Ma arises due to a common factor in the evaluation of the yAA form factors, which 
drops out of the ratios. One can speculate that the “quenching” effects drop out, at least 
partially, from the ratios as well. 

In Fig. HU we also show the mjr-dependence of the yAA transition ratios, with the the¬ 
oretical uncertainty estimated according to Eq. for the case LO= 0, and with the 

average taken over the range of from 0 to 0.15 GeV^. The dependence obtained here 
from yEFT clearly shows that the lattice results for Rsm niay in fact be consistent with 
experiment. 


VII. CONCLUSION 

Let us briefly go over the main points and results presented in this paper, which is the 
hrst one in a series devoted to yEFT in the A(1232)-resonance region. 

(i) We develop an extension of chiral perturbation theory to the A(1232)-resonance energy 
region, based on the ^-expansion of Ref. [1^ . In this yEFT framework the expansion 
is done in the small parameter 6 equal to the excitation energy of the resonance over 
the chiral symmetry-breaking scale. The other low-energy scale of the theory, the pion 
mass, counts as 5^, which is a crucial point for an adequate counting of the A-resonance 
contributions in both the low-energy and the resonance energy regions. 

(ii) This framework has been applied here to the process of pion electroproduction. This 
is a hrst yEFT study of this reaction in the A(1232)-resonance region. We have 
performed a complete calculation of this process in the resonance region up to, and 
including, next-to-leading order in the 5-expansion. The power counting in 5 has only 
been used to establish which graphs contribute at the leading and next-to-leading or¬ 
der, no actual expansion of the diagrams themselves is being done. Therefore, some 
higher-order in 5 effects, required by relativity and analyticity, are automatically in¬ 
cluded. Such effects are known to improve the convergence and extend the region of 
applicability of yEFT calculations. 
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(iii) Our NLO calculation of pion electroproduction satisfies gauge and chiral symmetries 
perturbatively, and Lorentz-covariance, analyticity, unitarity (Watson’s theorem) ex¬ 
actly. 

(iv) The chiral-loop contributions to the yiVA transition have been evaluated using two 
independent techniques: the Feynman parameter method and the sideways dispersion 
relations. Both methods yield the same result. 

(v) The only free parameters entering at this order are the 7 AA couplings qm, dE, gc 
characterizing the Ml, E2, C2 transitions, respectively. By comparing our NLO 
results with the standard multipole solutions (MAID and SAID) for the photopro¬ 
duction multipoles we have extracted gM = 2.9 and gE = —1.0, corresponding to 
Rem = —2.3 %. The NLO yEFT result was also found to give a good description of the 
energy-dependence of most non-resonant s, p and d-wave photoproduction multipoles 
in a 100 MeV window around the A-resonance position. From the pion electroproduc¬ 
tion cross-section we have extracted gc = —2.36, which yields Rsm — —7 % near 

= 0.1 GeV^. In overall, the NLO results are consistent with the experimental data 
of the recent high-precision measurements at MAMI and BATES. 

(vi) The yEFT framework plays a dual role in that it allows for an extraction of resonance 
parameters from observables and predicts their pion-mass dependence. In this way it 
may provide a crucial connection of present lattice QCD results (obtained at unphysical 
values of to the experiment. We have shown here that the opening of the A — eN 
decay channel at tUtt = Ma — Mn induces a pronounced non-analytic behavior of the 
Rem and Rsm ratios. While the linearly-extrapolated lattice QCD results for Rsm 
are in disagreement with experimental data, the yEFT prediction of the non-analytic 
dependencies suggests that these results are in fact consistent with experiment. 

(vii) The present calculation is systematically improvable. We have indicated what are 
the next-next-to-leading order effects, however, at present we could only estimate the 
theoretical uncertainty of our calculations due to such effects. We have dehned and 
provided a corresponding error band on our NLO results. An actual calculation of 
N^LO effects is a worthwhile topic for a future work. 


As high-precision data for low-Q^ pion electroproduction in the A-resonance region be¬ 
come available from BATES, MAMI and JLab, and the next-generation lattice calculations 
of the 7 AA transition are on the way 0 , the yEFT presented here makes a promise to be 
the theoretical framework to examine and connect these results. 


APPENDIX A: FURTHER EVALUATION OF THE FEYNMAN-PARAMETER 
INTEGRALS 

To perform the integrals in Eqs. (I46all46c|l . we hrst note that can be written as 
[{x-(5 + y(3.y)^-V^], with 

Riv) = [iP - yPif - - Q^y{^ -y) + , (Ai) 
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and then use the following elementary integrals 


dx ln[(a; — a)^ — 6^] = —2x + {x — a) ln[(a; — a)^ — b‘^] + b In 


X — a + b 
X — a — b 


+ C, 


dx 


dx 


dx ■ 


1 , X -a + b ^ 
= — In-r + C, 


{x — aY — 2b X — a — b 


X 


{x — a)^ — 6^ 


= ^ ln[(a; — a)^ — 6^] + 


a , X — a + b ^ 

Tb'''7^^b + ^’ 


(A2) 


X 


{x — a)2 — 6^ 
to perform one integration: 


= X + a ln[(a; — a)^ — 6^] + 


+ 6^ , X — a + b 


2b 


In 


X — a — b 


+ C, 


g‘f^ = -Cna I dyy I dx InM'i = -Cna I dyy {-2(1 - y) 


1 i-y 1 

r r r 

<2 


0 


+ (a - ya^) lnp(l - yf + yy^] + (/3 - yjd^) \n[y? + Q‘^y{l - y)] 
(1 - l/)(/5 - yP'y - yQ"^ + V)- + ie 


V In 


[{l-y){j3- yj3^ - yQ"^ -V)- y^ + ie 
1 ^-y 


Qe = Cna j dyy J dx (inAfg - 2x [-x + {1 - y) (1 + r)] = -Cna j dyy 

0 0 0 
X {{p-yp^- - y )) lnp(l - yf + yy^] + (a - ya^ + r(l - y)) \n[y^ + Q‘^y{l - y)] 

(1 -y){P- yP'y - yQ‘^ + v)- y^ + ie 


+ 


P-yPi 

V 


{a — ya^ + r(l — y)) In 


[{l-y){p- yp^ - yQ"^ -V)- y^ + ie 


(A3) 


^-y 


g'C = -Cna j dyy{2y - 1) j dx [-x + (1 - y) (1 + r)] M 
0 0 

rpl-yf + yy^ 


-2 

e 


= -Cna / dyy{2y- 1)1 


In 


-p(«-|/«7 + ’^(l-l/)) 111 


- y) 

(1 -y){.P- yP'i - yQ'^ + v) - y'^ + ie 


L(1 -y){p- yp^ - yQ"^ -V)- y^ + ie 
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Furthermore, define 


'In 

Jn 

ICn 


j dyy^ \n[{r / yf {I - yf + y], 

0 

1 

J dyy"^ ln[l + {Q/yfy{l - y)], 
0 


1 

j dyy^V In 
0 


(1 -y){(d- yp^ - yQ^ + V)- 
(1 - y)iP-yfd'r -yQ"^ 


[, y\ {l-y){P-yP,-yQ^ + V)-y^ 
J ^(l-y)(p-yP^-yQ2-V)-y^^ 


(A4a) 

(A4b) 

(A4c) 

(A4d) 


Then, 

g^jJ = Cj\fA { 5(1 — In/i) — crXi + 0 :^X 2 — /3j7i + — /Ci} , (^5) 

= -Cna {-^\ny- {P - r)Xi + (/3^ - r)X 2 - (a + r)Ji + (a^ + r)j 2 

+ /3(q; + r )Xi — [/3(ci^ + r) + j3^{oi + r )]X2 + + r )X3} , (-^6) 

9c^ ~ ^C'ata { 2 X 2 — Xi — 2J72 + iXi + (ci + ^")(Xi — 2 X 2 ) ~ (CI 7 + ^")(X 2 — 2 X 3 )} , (A7) 


Xi = 


Xo = 


8 - 6Wr)^ + (^Irf ^ ^ 

V4(r/A.)'" - 1 

+ \{y/rf + [1 - 2{y/rf + \{y/rf] \n{r/y) , 
^-^i{ylrf + Uylrf-Uy/rf 


VWIW 


■ arctan 


Ji 

J2 


= -1 + 


f + I- U^^/rY + [f - 3(/i/r)2 + 2(/i/r)^] \n{r/y) , 
v/ 4 /i 2 + g2 ^4^2 + g2 + g 


2g 


v/4/i2 + g2 _ g ’ 


^ _ 2 ^ ^ + -Q + '^) A/4/i2 + g 2 + g 

18 3g2+ 3v/V + g2 "" ^Ay^ + - Q ' 


(A8a) 


(A8b) 

(A8c) 

(A8d) 


APPENDIX B: ABSORPTIVE PARTS OF THE VERTEX CORRECTIONS 

In the sideways dispersion relation method, we evaluate the ttN loop contribution to the 
'jNA form factors using a dispersion relation in the invariant mass of the A. To present the 
results of this technique, we first introduce a number of kinematical quantities and notations. 
The pion momentum in the rest frame of the decaying A with invariant mass W = yX is 
given by Ip^I = hF A, with 

A = ix/[l - (r+ /f)2] [1 - (r-/f)2], 

= [{l-fy-~y^], (Bl) 
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and fi = r = Mn/W. We also introduce ta = Mj\/W. In the on-shell case, 

W = Ma, we of course have f = r, p. = fi, and A = A. 

Furthermore, the energies of the pion and nucleon (-Eat) in the rest frame of the 
decaying A can be expressed as = W P and = Wa with the fractions given by: 

P = ^ [l - , a = ^ (l + - p"^) = 1 - p. (B2) 

and the Lorentz-invariant pion (nucleon) velocities, denoted by {vn) respectively are given 
by: 


= 


IPtt 

E^ 


A 

f 


Vn = 


IPtvI 
E 


N 


A 

a 


(B3) 


where we used |p 7 v| = IPttI in the A rest frame. 

We denote the virtual photon four-momentum hy q = {q^,<^)- For spacelike virtual 
photons, the virtuality is usually denoted by — gg > 0. The virtual photon three- 

momentum in the A rest frame can then be expressed as |q| = W R, with 


K = 


(I_y2)2 + 2g2(i + y2)+g4^ 


(l + h)2 + Q2j ^(l_y)2 + Q2j ^ 

where js. It is also useful to introduce the fraction of the photon energy: 


P^ = 


V -q 


= h^-r^-Q' 


satisfying Finally, we dehne 


2 kX 


= 


Wn = 


2 IpttI |q| __ 

g2 + 2E^gO “ Q2 +2/3/3.’ 


2 IpjvI |q| 


2 kX 


(B4) 

(B5) 

(B6) 

(B7) 


g2 2 E]\f g2 -j- 2 (5 /3g, 

In the real photon limit 0), one obtains that w.^ —> and wn —> v^. 

The calculation of the absorptive part of these ttA loop contribution, entering the side¬ 
ways dispersion relations, involves several phase space integrals, of which the simplest is of 
the form : 


[(2vr)5(p2-m2)0(pO)] [(27r) (5((p - - M^) 0(p'° - p°)] = — 


dvr 


(B8) 


The evaluation of the absorptive part corresponding with Fig. El (a), hrstly involves the 
following phase space integrals with different powers of pion momentum in the numerator : 


/ [(27r)5(p^-m2)0(pO)] [(27r)(5((p'-p^)2-M^)0(p'°-p°)] p" 

= ai,^p'“, (B9) 


/ [(27r)^(p^-m2)0(pO)] [{2tt)6{{p'- p^p 

= a2,^p'^g^^ + a3,.p'>'^ 


Ml) 0(p'° -p°)] p> 


A 

TT 


(BIO) 
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with 



A ;; 

,7r 



^ ^ ( 

^2 ,7r 

4 vr y 


1 r 

,7r 

— — A 

4 vr 


3 




(Bll) 

(B12) 

(B13) 


Furthermore, the calculation of the absorptive part corresponding with Fig. El (a), which 
has one pion which is off-shell, also involves the following integrals : 

f 


(2i)‘ 


2p‘ 


d'^PiT 

(27r)‘ 


2p 

d‘^p^ 

W) 


[{2'k)5{pI 


,/2 


[{2ti)6{pI 



+ &3, 

[(27r)^(p2 



T - 2p.T ■ q 


(B14) 


[{2'k)5{pI - ml) 0(p°)] [(27r) 5{{p' - p^)^ - M^) 0(p'° - p°)] 




T - 2p.T ■ q 


(B16) 


pIpI 


If a/3, / (q°p'^ + q^p'’^\ , C2,^ , C3,^ „ p\ 

+ j+^P P' + TT?/}. 


g 2 - ■ q 

(B16) 


/ p7r)S{pl-ml)e{p°)] [(2^)i((p'-p,)=-Mpe(p'»-p;)] ^ 

-h d\ ^ rg"p'+ +gW^l 

' q- pf ^ j 

+ d2,. 

+ 4 , ^ 

’ q ■ p ^ j 

. ,7r /a / 3 /■y , ,7r a /3 a 1 

+ P^P^ + —^q q^q^ > ■ 

ip' ^ ip' ^ I 


(B17) 
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The coefficients appearing in the above equations are given by : 


T± log 

An 2k \1 — Wt, ) 


1 ^ f o 1 . 


N k 


1 + Wt. 
1 - ta-7, 


1 + tc„ 

l-W-n 


(B18) 

(B19) 


(B20) 


—-|-T + log('i±^ 

An Ak \ Wn V 1 — tc„ 


-1 + 


An Ak^ 1 ( 


1 (l + w^ 

k‘°® 


i_J_JA h+i^h 

AnAk^ 1 \ 2 k^ 


2 ( 3 ( 3 ^ + ) + 2 Q 




g^(2/3/3, + g^) ^ (B22) 


+ ) — 4Q^ 


^logfi±^ 

2k ^\l-w. 


i _^ , A + 

47r 1 ^ \ 1 — 


2 AM 1 - — + 4/^2 + g 


-1 + 


(B21) 


(B23) 


(B24) 


1 I 4,^5 QH2P - n) , Q\ /'l + KI,^,„,5 
4^53 - 3 ''^'’’- n -+ T ‘“8 rv^r (2/3 - /S.) 


-1 + 




+ 2Q^ + 2(3(3. 


2 k V 1 - w. 


A" (L(2fi - A) h - 4V k" fl - A 


+/3(2/3 — (3.y)(2(3(3.y + Q‘ 


, (B25) 


(B26) 
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^2 .7r 


1 r2 


77-? + 77— log 


477 2^2 1 3 w‘i 2m 


d' = — _ — 

2’" 477 8^4 
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2A^ -77 + 


1 + 
l-w^ 
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(B27) 


4k"/3/3^ 1 


wt 


(3,Q\2(3-(3,)[3 


wt 


d^,-K = ~3 C?1 ^TT — 4 — 2/3^((i2,7r + C^2 , 7 r) +^2,J 

2 1 

<^4,7r = ^('^2,7r + C?2,7r) + 77 03,7r- 

g"= / 


(B28) 

(B29) 

(B30) 


The calculation of the absorptive part corresponding with Fig. El (b), which has one 
nucleon which is off-shell, involves the integrals : 


d'^p 


'N 


(277)4 

ai,. 

d^^PN 


[(277) 6{p% - Ml) 0(p^)] [(277) 5((p' - pnY - ml) 0(p'° - p^)] pi 


= apvrp'", 


(B31) 


[(277) 5(p^ - M^) 0(p^)] [{2'k)5{{p -pnY -ml) -pD] pi pi 


J (277)4 

= a2,NP'^g^^ + as^NP'^^p'l 

[(277) (5(p^ - M^) 0(p^)] [{2Tr)6{{p' -pn)^ -mDQip''^- pD] - ^ 


(B32) 


2p‘ 


/2 


hi 


N, 


g 2 - 2pN ■ q 
(B33) 


and analogous integrals as in Eqs. (IB15IIB17F where p^ ^ Pn- 
The coefficients appearing in the above integrals are given by : 

A _ 


Ol ,7V 

4ir“’ 




A3 

0^2,N 

4 77 Y 

3 

0^3,N 

= ^A 

4 77 



4A2 


and 


i. 1 1 1 fl + WN^ 

bi,N = ^7777 log -- , 

477 2 K yl- — '^N J 

h 

4774^2 


, 1 1 


/I + WatA 1 


4A/3,+ Q4 i + y2 + Q2 + 2(r~2-p2) log^^ 

L J 2 K — wn ) J 


-2A + 


+ g^ 77-log 


/I +Wn\ \ 


2k Vl -m^y J 


(B34) 

(B35) 

(B36) 

(B37) 
, (B38) 
(B39) 
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Cl ,N — “{ — 2 / 3 ^ai ^AT + 2 (r — Q \ j bi pf + Q f3^b2 + Q (^ + P'y) b^ 
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i2A%„ 
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C2,Ar — —2 Cl ,Ar — C3,Ar + 6l ,Ar + ^3,Ar, 
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C3,N — 




} (B40) 

(B41) 

(B 42 ) 

(B 43 ) 


di,N — 

+ 


2^2 

2/3^ 


—(2 a2 ,Ar + as ,7v) + — Q^) — 


^2 ,N 


2R^ -Q^ 


^1 ,N ^ id'y (C 2 ,Af + 2 Cl ,Ar) ^ , 


(B 44 ) 


5 , , (1 — r^) 


^\,N — C?1,7V + — &2,Af — 
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51 

2^2 


4^2 


a3,Af 


5 (3i ,Ar + - ^3 ,Ar — 62 ,Ar + -/37 (c'j^ + C 2 ,Ar) 


, _ 1 3 , 1 1 

“2,AT — /37 Cti^Ar + - Ci,Ar — ^ a2,Af, 

*32,AT = 2 g 2 ,Af + 2 <3i ,Ar) 3" 4 *^1 ,A^’ 
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2 1 

<34,Ar = ('32,Af + <32,Ar) + o C3,Af- 

g2 ^ 


(B 45 ) 
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(B 47 ) 
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The imaginary parts of g\^ g 2 , and g^ corresponding with the diagram of Fig. El (a) are 
thus given by: 

T- + dAh-A Mv 
Im^l,a = 




4J2 {d2,n} , 

gA Ha Mai ^ J (r a + f) 

W 




4/2 ^ 

gA hA Mn J_ J (r a + r) 
4/2 fT\ 2 


/ r, ^1 J 1 I ^ ^ ^ 

(2 <^2 ,7r ~ C;^ 




(B 50 ) 


(&3 , 7 r - 3 C3 ,^ + 2 (i4 (4 d '2 - C4 
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The imaginary parts of gi, g 2 , and g^ corresponding with the diagram of Fig. El (b) are 
likewise given by: 


Imp 6 = 
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Impfe = 


gA hA Mn 
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gA hA Mm 
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FIG. 1; Nucleon self-energy to one loop. 



FIG. 2: Diagrams for the eN —> eirN reaction at NLO in the ^-expansion, considered in this work. 
Double lines represent the A propagators. 


fi 

q 

p p' 

FIG. 3: 7 *iVA vertex. The four-momenta of the nucleon (A) and of the photon are given by p {p') 
and q = p' — p respectively. The four-vector index of the spin 3/2 field is given by a, and p is the 
four-vector index of the photon field. 
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FIG. 4: The higher-order 'jNA vertex correction required by the electromagnetic current conser¬ 
vation. 



FIG. 5: Absorptive part of the ypA"*" vertex. Diagram (a) : vr^n loop where the photon couples to 
the Tr'*'; diagram (b) : loop where the photon couples to the charge of the proton. The vertical 

dotted lines indicate that the ttN intermediate state is taken on-shell. 





FIG. 6: Examples of N^LO contributions to the eN —> eirN reaction neglected in this work. Open 
circles denote the e.m. coupling to the anomalous magnetic moment of the nucleon. 
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f3/2) f3/2) 

FIG. 7: (Color online) Multipole amplitudes (left panel) and (right panel) for pion 

photoproduction as function of the invariant mass W of the t:N system. Green dashed curves: 
A contribution without the yAA-vertex corrections, [i.e., Fig. Oja) without Fig. He, f)]. Blue 
dotted curves: adding the Born contributions. Fig. Hb), to the dashed curves. Black solid curves: 
complete NLO calculation, includes all graphs from Fig. ^ In all curves the low-energy parameters 
are chosen as : qm = 2.9, gE = —1.0. The data point are from the SAID analysis (FA04K) 

(red circles), and from the MAID 2003 analysis [ll| (blue squares). 
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FIG. 8: (Color online) Non-resonant multipole amplitudes (in units 10“^/m7r) for pion photopro¬ 
duction as function of the invariant mass W of the ttN system. The solid curves result from our 
NLO calculation. The data points are from the SAID analysis (FA04K) (red circles), and from 
the MAID 2003 analysis (blue squares). 
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FIG. 9: (Color online) The NLO results for the 0,^ dependence of the 
at W = 1.232 GeV and = 0.127 GeV^. The theoretical error bands 
Data points are from BATES experiments iQ. 


7 *p —> TT^p cross sections 
are described in the text. 
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FIG. 10: (Color online) Non-resonant multipole amplitudes (in units 10“^/m7r) for pion electro¬ 
production at the resonance position {W = 1.232 GeV) as function of Q^. The curves are the 
results of our NLO calculation. The data points are from the SAID analysis (FA04K) (red 
circles), and from the MAID 2003 analysis |lll| (blue squares). 
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FIG. 11: (Color online) dependence of the NLO results (solid curves) for Rem (upper panel) 
and Rsm (lower panel). The blue dashed curves represent a phenomenological estimate of N^LO 


effects by including Q^-dependence in 


from MAMI for Rem 111, and Rsm 


4C 


a ccording to Eq. (PH) . The blue circles are data pomts 
. The green squares are data points from BATES 
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FIG. 12: (Color online) Pion-mass dependence of the nucleon and A(1232) masses. The curves are 
two-parameter expressions for the ttN loop contributions to and Ma according to Eqs. (HU) 
and (0^ respectively, using = 0.883 GeV, cim = —0.87 GeV“^, and 
—0.40 GeV~^ respectively. The red squares are lattice results from the MILC 
The stars represent the physical mass values. 


= 1.20 GeV, ci^ 
Collaboration I 43 I. 
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FIG. 13: (Color online) The m^-dependence of the real part of the Jones-Scadron 'yNA form 
factor for = 0 and = 0.127 GeV^. The solid (dashed) curves are the NLO results for 
= 0.127 GeV^ (Q^ = 0) respectively, including the dependence of M^r and Ma. The green 
dotted curve is the corresponding result for = 0.127 GeV^ where the mjr draendence of Mn 
and Ma is not included. The blue circle for = 0 is a data point from MAMI []| , and the green 
square for = 0.127 GeV^ is a data point from BATES [^. The three filled black diamonds at 
larger m-,^ are lattice calculations for values of 0.125, 0.137, and 0.144 GeV^ respectively, 
whereas the open diamond near m-,^ ~ 0 represents their extrapolation assuming linear dependence 
in m^. 
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FIG. 14: (Color online) dependence of the NLO results at = 0.1 GeV^ for Rem (upper 
panel) and Rsm (lower panel). The blue circle is a data point from MAMI , the green squares 
are data points from BATES j^. The three filled black diamonds at larger m,r are lattice calcula¬ 
tions [ 2 ^, whereas the open diamond near ~ 0 represents their extrapolation assuming linear 
dependence in m^. Red solid curves: NLO result when accounting for the dependence in M^v 
and Ma; green dashed curves: NLO result of Ref. 0 , where the m,r-dependence of Mjv and Ma 
was not accounted for. 
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